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ABSTRACT 

The distance-redshift relation plays an important role in cosmology. In the standard 
approach to cosmology it is assumed that this relation is the same as in the homo- 
geneous universe. As the real universe is not homogeneous there are several methods 
to calculate the correction. The weak lensing approximation and the Dyer-Roeder 
relation are one of them. This paper establishes a link between these two approxima- 
tions. It is shown that if the universe is homogeneous with only small, vanishing after 
averaging, density fluctuations along the line of sight, then the distance correction is 
negligible. It is also shown that a vanishing 3D average of density fluctuations does 
not imply that the mean of density fluctuations along the line of sight is zero. In this 
case, even within the linear approximation, the distance correction is not negligible. 
The modified version of the Dyer-Roeder relation is presented and it is shown that 
this modified relation is consistent with the correction obtained within the weak lens- 
ing approximation. The correction to the distance for a source at z ~ 2 is of order 
of a few percent. Thus, with an increasing precision of cosmological observations an 
accurate estimation of the distance is essential. Otherwise errors due to miscalculation 
the distance can become a major source of systematics. 
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1 INTRODUCTION 



gravitational lensing: weak - cosmology: theory - large-scale structure 



against cosmological observations by Santos & Lima (2006), 
Santos, Cunha & Lima (2008), and Yu et al. (2010). 



The distance-redshift relation plays an important role in cos- 
mology. In fact almost all cosmological observations depend, 
either explicitly or implicitly, on this relations. However, in 
general case, to calculate the distance, knowledge of matter 
distribution (as well as its evolution) between the observer 
and the source is needed. Therefore, an approximate rela- 
tions are of great use to astronomers. The simplest approx- 
imation assumes homogeneity of the universe and is based 
on the Friedmann model. However, as the real universe is 
not homogeneous a more elaborate relations are needed. 
Zel'dovich (1964) proposed an approximation which takes 
into account that light propagates through emptier rather 
than denser regions of the universe. This approximation 
is now known as the Dyer-Roeder approximation (Dyer & 
Roeder 1972, 1973). It also assumes homogeneity but allows 
for a different density than in the background model. The 
density difference is modelled by a constant, the so-called 
smoothness parameter a. As the Dyer-Roeder equation is a 
differential equation an approximate analytic solution was 
presented and discussed by Demianski et al. (2003). Gener- 
alisation to a(z) was first suggested by Linder (1988) and 
the effect of the change of the expansion rate was discussed 
by Mattsson (2010). The Dyer-Roeder relation was tested 



However, it has been argued that the Dyer-Roeder re- 
lation may not properly describe the effect of matter clus- 
tering (Rasanen 2009; Ellis 2009) and therefore may not be 
an appropriate approximation for the distance-redshift rela- 
tion in the real universe. Also, if the smoothness parameter 
changes with redshift the distance formula depends on sev- 
eral free-parameters whose interpretation can be ambiguous. 
Therefore, although cosmologists are aware that not taking 
into account inhomogeneities introduce additional system- 
atics, the Dyer-Roeder relation is not widely used. Instead, 
several alternatives based either on the linear perturbative 
scheme (see Bonvin, Durrer & Gasparini 2006 and refer- 
ences therein) or non-linear models (see Bolejko et al. 2009 
and references therein) have been developed. 

This paper explores the connection between the Dyer- 
Roeder approximation and the weak lensing approximation. 
It shows that the Dyer-Roeder relation can be modified so 
that it is consistent with the lensing approximation. The 
modified version contains only 1 free parameter of clear in- 
terpretation, i.e. the mean of density fluctuations along the 
line of sight {8) id- 

The structure of this paper is as follows: Sections[2]and[3] 
present the Dyer-Roeder and weak lensing approximations 
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respectively, Sec. [4] presents the comparison of these two 
methods, and Sec. [5] discusses the results. 



2 THE DYER ROEDER APPROXIMATION 

The angular diameter distance Da is given by the following 
relation (Sachs 1961) 



A 2 D A 
da 2 



(1) 



where a is the shear of the light bundle, k a is a vec- 
tor tangent to the light ray, R a p is the Ricci tensor, and 
R a pk a k^ = KT a pk a k^ (where T a p is the energy-momentum 
tensor). In the comoving and synchronous coordinates, for 
pressure- less matter, T a fjk a k^ — Kpk°k°. The Dyer-Roeder 
approach assumes homogeneity (a — = Sp) but takes into 
account that light propagates through vacuum. Therefore, 
Q, m which photons 'feel' is different than true fi m . This is 
modelled by a constant parameter a (of value between and 
1) the multiplies fi m . In this case {TJ reduces to 



d D a 

dz 2 



- + 



/ 1 dH 
\H dz 



1 + z 



dD A , 2aQ, m Hl 



dz 



- + 



2H 2 



(l+z)D A = 0, (2) 



where H{z) = H ^/Q, m {l + z) 3 + Q k (l + z) 2 + fi A . The 
initial conditions needed to solve ([2j are: Da = and 
dD A /dz = 1/Ho. 

Generalization to a(z) was explored by Linder (1988) 
who suggested several algebraic forms like a(z) = aa + aiz 
or a(z) = q» + 02(1 + z) 1 . Santos & Lima (2006) proposed 
a(z) = /3 (l + z) 37 /[l + A)(l + 2) 37 ]- This paper studies the 
following form of a (2) 



x(z) = 1 + T>(z)(8)i 



(3) 



where (5) id is the mean of present-day density fluctuations 
along the line of sight, and T>(z) describes its evolution. Be- 
low it is shown that if V = (1 + z) -5 / 4 , then the Dyer- 
Roeder equation gives results consistent with the results ob- 
tain under the assumption of the lensing approximation. 



LINEAR PERTURBATIONS AND THE 
LENSING APPROXIMATION 



Writing the distance as 
D A {z) = D a (1 + 5d), 



(4) 



where Da the distance in the homogeneous universe, one can 
derive formula for 8d using the linear perturbative scheme. 
The most general form was presented and discussed by Pyne 
& Birkinshaw (2004), Bonvin et al. (2006), Hui & Greene 
(2006), and Enqvist, Mattsson & Rigopoulos (2009). Ex- 
cluding the contribution from the motion of the observer 
and source, and taking the leading term, 8d reduces to 



Sn = - 



Xe 



Xe 



x -xv 2 0(x), 



(5) 



where \ i s the comoving coordinate dx = dz/H(z), <j> is 
the gravitational potential which can be related to density 
perturbations p8 via the Poisson equation X7 2 (p = ^£-a 2 p5. 
Equation J5} is equivalent to the convergence in the lensing 



approximation and is known as the Born approximation. As 
seen voids (<5 < 0) increase the distance while regions of 
S > decrease it. 

In order to solve ([5]) one needs to know V 2 <f> along the 
line of sight. Using the Poisson equation, the gravitational 
potential is related to the density fluctuations. Thus when 
calculating the variance, one can Fourier transform S and use 
the matter power spectrum instead (Munshi & Jain 2000, 
2001). However, in this paper we are not interested in the 
variance or higher oder moments, just in the distance itself. 
Therefore, to solve (f5} actual density fluctuations along the 
line of sight are needed. 

The preset-day fluctuations are non-linear and there- 
fore no longer Gaussian. Only when density fluctuations are 
in the linear regime, their probability distribution function 
(PDF) can be approximated by the Gaussian PDF. Once 
density fluctuations are in the non-linear regime they are no 
longer Gaussian (there is no symmetry as — 1 ^ S < 00). 
However, it has been shown that in the non-linear regime 
density fluctuations can be approximated by the one-point 
log-normal PDF (Kayo, Taruya & Suto 2001; Lahav & Suto 
2004) 



P(S) = 



v/2^ 



exp 



(ln(l + 5) + ^ ; /2) 2 



where 



Ml + and o R 



1<l 



1 

2^2 



1 



1 + 5' 



(6) 



dkV(k)W 2 (kR)k 2 ,(7) 



where V(k) is the matter power spectrum. 

Two methods for generating density fluctuations along 
the line of sight are considered. Both have a log-normal PDF 
of 8 (3D). The first method ensures that (<5)izj ~ (the mean 
of the present-day fluctuations along the line of sight). The 
second method, on the other hand, allows for {5} m 7^ 0. 



4 RESULTS 

4.1 Vanishing mean of the density fluctuations 
along the line of sight 

First let us focus on generating the density fluctuations using 
directly the log-normal PDF. Using ([6]) 5 smoothed in a 
sphere of radius R can be generated. Density fluctuations in 
this model are schematically presented in Fig. [TJl — when the 
light ray exits one sphere it enters another one of different 
R and S (for details see Appendix 1X1. 

The distance correction 5d calculated from ([SJ is pre- 
sented in Fig. [2] One point about Fig. [2] needs to be empha- 
sized. As follows from ([5| the distance correction depends 
on the position of the source, for example 5d{z = 0.5) for a 
source at z» — 2 is different than 5d(z = 0.5) for a source at 
z* = 0.5 even if density fluctuations along the line of sight 
are the same. Thus, the distance correction presented in Fig. 
[2] is the distance correction for the source at given z, i.e. 
8n{z = z- t ). Fig. [2] presents the mean and variance: at each 
z (a discrete number of z with an interval A2 = 0.01 was 
considered) based on 100,000 runs (each with a different dis- 
tribution of density fluctuations along the line of sight) the 
mean and variance were calculated. The mean and variance 
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Figure 1. Density inhomogeneities along the line of sight. Panel 
a): density fluctuations generated from the log- normal PDF (0 - 
when the light ray exits one structure, the next one is generated. 
Panel b): compensated density fluctuations - when the light ray 
exits a structure, another one is generated together with an initial 
angle at which the ray enters the structure. When 7 is large then 
the light enters and exits the structure at large angles. In this 
case structures can overlap. Panel c): light propagation in the 
real universe. 



0.006 r 
0.004 - 



0.002 - 
- 
-0.002 - 




-0.004 - 



-1 -0.5 0.5 1 1.5 2 

5 

Figure 3. PDF of <5 smoothed in a sphere of radius 4 Mpc for a 
model discussed in Sec. I4.2I (solid line) and the log-normal PDF 
(|6j of 8 smoothed also in a sphere of radius 4 Mpc (dashed line). 
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Figure 4. Mean density fluctuations along the line of sight (O 
for a model discussed in Sec. I4.2I The upper most curve presents 
a model with 7 generated from a uniform distribution between 
and 0.5 it. The second curve presents a model with 7 generated 
from the Gaussian distribution with the mean equal to zero and 
the standard deviation equal to 0.25-7T. In the third model 7 is 
generated from a uniform distribution between and 0.257T. The 
bottom curve presents a model with 7 = 0. 
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Figure 2. Distance correction 8d within the weak lensing ap- 
proximation for a model with a negligible mean of density fluctu- 
ations along the line of sight (Sec. l4.iT l. Dash-dotted line presents 
80 obtained using the Dyer-Roeder relation with the smooth- 
ness parameter a = 0.99 (a = 1 is equivalent to the standard 
Friedmann relation). 

are represented by solid lines. For comparison, the Dyer- 
Roeder approximation with a = 0.99 is also presented (when 
q = 1 the Dyer-Roeder approximation reduces to the stan- 
dard Friedmann relation for the distance). As seen 8d is of 
negligible amplitude ~ 10 -3 . Thus, for this configuration, 
and under considered assumptions, there is no need to take 
into account inhomogeneities. Also, another important fact 
is that Sd in the weak lensing approximation is negative, 



while the Dyer-Roeder approximation implies positive 5d- 
This is a consequence of positive density fluctuations along 
the line of sight. For example, for a source at z — 3 (based 
on 100,000 runs) the mean of the present-day density fluc- 
tuations along the line of sight defined as 

r 

(5) 1D (r) = ± J df5(t ,f), (8) 


is (5) id — 6 x 10 -3 and the standard deviation 2 x 10 -2 , 
which implies 3d = —5 x 10~ 4 ± 3 x 10~ 3 . Thus, as pointed 
out by Rasanen (2009) in order to properly describe the ef- 
fect of clustering, the smoothness parameter should be larger 
than 1, not less than 1 as in the Dyer-Roeder approxima- 
tion. However, in this case the difference as well as the total 
correction are negligible. However, even if the mean is negli- 
gible, the variance can provide additional information - see 
Linder (2008) for a discussion. 
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4.2 Non-zero mean of the density fluctuations 
along the line of sight 

In the previous section, random matter fluctuations along 
the line of sight were considered. In such a case {5)3d ~ 
implies (<5)id = 0, and the distance correction is negligible. 
However, present-day matter fluctuations in the Universe 
are not purely random but are organized - matter in the 
Universe forms the cosmic web. In this case {8)30 = does 
not necessarily imply (5) id = 0. As seen in Fig[TJ; the cos- 
mic web contains large voids with fairly compact filaments. 
Therefore, photons spend more time in voids than in over- 
dense regions. In order to model this phenomenon let us 
consider a universe that consists of voids surround by fila- 
ments, each structure compensated. Thus, by the construc- 
tion a 3D average of density fluctuations is zero (if averaged 
over sufficiently large scale). Moreover, if the parameters of 
the system are adjusted (see Appendix [B] for details) then 
the PDF of density fluctuations can be almost log-normal. 
This feature is presented in Fig. [3] which shows the PDF 
of density fluctuations smoothed within a sphere of radius 4 
Mpc. As seen the PDF is similar to log-normal PDF, apart 
from high 5 where PDF is of lower amplitude - but see Jain, 
Seljak & White (2000) where the PDF of density fluctu- 
ations smoothed on scales of 3 h' 1 Mpc has also a lower 
amplitude for high 8 than the log-normal distribution. 

When the light ray exits one compensated structure 
before it enters another one, not only parameters of the 
next structure are generated, but also an angle 7, at which 
the light ray enters another structure. Therefore, 4 different 
methods are going to be considered 

(i) 7 = 0, the light ray always enters another structure at 
7 = 0, thus it passes through the centre of the structure. 

(ii) 7 is randomly generated from a uniform distribution 
between and 0.5 tt. However, as shown in Fig. [Tp in this 
case structures can overlap. 

(iii) to reduce overlapping, 7 is generated from the Gaus- 
sian distribution with the mean and a = 0.257T (if | — )/ 1 > 
0.5-7T then I7I — > n — 

(iv) 7 is randomly generated from a uniform distribution 
between and 0.25 n. 

The distance correction Sd (the mean and variance) is 
presented in Fig. [5] As seen if 7 is randomly generated from 
a uniform distribution between and 0.5 n then the mean 
of density fluctuations along the line of sight is small and 
Sd is almost negligible (cf. Brouzakis, Tetradis & Tzavara 
2008; Vanderveld, Flanagan & Wasserman 2008). 

However, as pointed out above, if 7 changes between 
and 0.57T then structures can overlap (see FigQ~p). To reduce 
the overlapping, 7 needs to by chosen from a smaller range 
of angles. In this case 5d increases. This is because the light 
rays propagate more likely through voids than filaments. 

In Fig. [S] the Dyer-Roeder relation is also presented. 
The dash-dotted line presents Sd obtained from the origi- 
nal Dyer-Roeder formula, the dashed line and dotted line 
prenent the evolving smoothness parameter ©. Intuitively, 
if {S)id were just an ordinary density perturbation, in the 
linear regime, it should evolve according to (|A1|| . Thus, the 
dotted line presents a model with a(z) — 1 + 5(2), where 8(z) 
is given by (1A1[) (with an initial condition S(z = 0) = (5)io). 

Following Mattsson (2010) let us consider the Dyer- 



Roeder relation with a perturbed expansion rate. From the 
continuity equation p + 3Hp = 0, the perturbation in the 
expansion rate are AH = —S/3 (where S follows from (|A1[I 
with an initial condition S(z = 0) = (8) id)- We use this 
relation to calculate the perturbed H(z) and AH/Az, and 
we then insert it to and using a(z) = 1 + S(z), where 
5{z) is given by (|A1|) . we get the distance correction that is 
presented using the dashed lines. 

As seen, neither a constant a, nor the evolving one (with 
8(z) given by (|A1[| ) . or with a pertubed expanion rate, are 
consistent with the weak lensing approximation (expect for 
the case where Sd is almost negligible). However, if a(z) = 
1 + (5)id/(1 + z) 1 , where 7 = 5/4 then the Dyer-Roeder 
approximation produces results comparable with the weak 
lensing approximation. As seen, the evolution of a is not 
as intuitively expected, i.e. it does not directly follow from 
(|A1[) . Here an empirical approach was employed, and it was 
found that if 7 ~ 5/4 then the modified, in this way, Dyer- 
Roeader relation leads to the agreement with the lensing 
approximation. 

An important result of the above analysis is that a non- 
zero mean of density fluctuations along the line of sight can 
modify the distance by a few percent. As cosmological ob- 
servations are now reaching the precision of a few percent, 
a proper handling of the distance is essential. Otherwise the 
errors in the estimation of the distance can become a ma- 
jor source of systematics. For example a proper handling of 
distance will be of great importance in the future analysis 
of BAO experiments (Bolejko 2010). 



5 CONCLUSIONS 

The distance-redshift relation plays a central role in cos- 
mology. If the homogeneous Friedmann model correctly de- 
scribes the evolution of the universe on large scales and, in 
addition, if density fluctuations along the line of sight van- 
ish after averaging, then within the linear approximation, 
the distance correction is negligible - it is sufficient to apply 
the Friedmann relation. However, if only the mean of den- 
sity fluctuations along the line of sight is not zero then, even 
within the linear approximation, the distance can change by 
several percent. 

In Sec. 14.21 it was shown that a vanishing 3D average of 
density fluctuations does not imply that the mean of den- 
sity fluctuations along the line of sight is zero. It is argued 
that in the real universe this may be the case. In the real 
universe voids occupy large regions while overdensities are 
more compact. Moreover, if light propagates for a long time 
through filaments then it is more likely to be absorbed or 
scattered. Thus, if a remote galaxy is observed then most 
likely its photons propagated through emptier rather then 
denser regions. In this case the weak lensing approximation 
produces results that are similar to results obtained from the 
modified version of the Dyer-Roeder equation. The modified 
relation is 

A 2 D A f±dH 2 \ AD A 
Az 2 \H dz l + z) Az 

l§n m (i +z) (i + T ^) DA = a, (9) 
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Figure 5. Distance correction So for a model discussed in Sec. I4.2I The weak lensing approximation J5} is presented by thin solid curves 
(mean + variance), the dash-dotted lines present the Dyer— Roeder approximation with a constant smoothness parameter a = 1 + (<5)i£>, 
the dotted lines represent models with ct(z) = 1 + S(z), where <5(z) is given by dAlll , the dashed lines present the distance correction for 
models with a perturbed expansion rate (see text for details), and the thick solid lines represent the modified Dyer-Roeder relations with 
a = 1 + (S)m /(l + z) 5 / 4 . Panel a): 7 generated from a uniform distribution between and 0.5 7T, (S)id ~ —0.04. Panel b): 7 generated 
from the Gaussian distribution with the mean and a = 0.257T, (&)\d ~ —0.19. Panel c): 7 generated from a uniform distribution 
between and 0.25 n (5) 1D fa -0.31. Panel d): 7 = 0, (8) 1D Ri -0.42. 



(with initial conditions Da ~ and ADa/Az — I /Ho). Thus 
apart from the background cosmological parameters, the 
only free parameter left is the mean of density fluctuations 
along the line of sight (8) m- Since ([9| is just an ordinary 
differential equation, it is as easy to implement it and solve 
numerically as the standard relation for the distance in the 
Friedmann model. The mean of density fluctuations along 
the line of sight (5) id could either be deducted from the 
galaxy redshift surveys or N-Body simulations. 

It should be noted that the analysis presented here was 
based on the linear approximation (the lensing approxima- 
tion), i.e. higher order corrections to the distance were not 
taken into account. It is likely that quadratic corrections 
may lead to non- negligible So even when (5) id ~ 0. In 
this case the presented above modified version of the Dyer- 
Roeder equation would not be consistent with the actual 
non- linear distance-redshift relation, as when (<5}id = it 
reduces to the Friedmannian relation.. 

The most important conclusion of this paper is that 
even within the idealized case of large scale homogeneous 
universe (with only Mpc-scale inhomogeneities) the distance 
can be different than in the homogeneous universe. The dif- 
ference can be of order of a few percent. Thus, with an in- 



creasing precision of cosmological observations an accurate 
estimation of the distance is essential. Otherwise errors due 
to miscalculation the distance can become a major source of 
systematics. 
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APPENDIX A: ALGORITHM OF 
CALCULATING 8 D FOR MODEL DISCUSSED 
IN SEC. 4.1 

The algorithm of calculating 5d along the line of sight is as 
follows 

(i) The radius of a structure R is randomly generated 
from a uniform distribution, from to 3 Mpc. 

(ii) Then an is calculated using ®. The primordial 
power spectrum was chosen in agreement with the WMAP7 
data (Komatsu et al. 2010): Afc" with n s = 0.969 and the 
amplitude A chosen so that erg = 0.803. The transfer func- 
tion was calculated according to Eisenstein & Hu (1998). 

(iii) An initial value of a density fluctuation So is gener- 
ated from the log-normal distribution ([6]). 

(iv) The evolution of S (at a fixed point) is calculated 
using the linear approximations (Peebles 1980) 

■f „ d i 4irG r , . , i 

<5 + 2-5 = — 5-p5- (Al) 
a c 

(v) Using the Poisson equation, V 2 is calculated and 
inserted into (J5j) which is solved from \i to x — X* + x(2-R)- 



(vi) Steps (i)-(v) are repeated so ([5]) is solved from x = 

to Xe- 

The cosmological parameters are the same as derived 
from the 7-year WMAP data: H = 71.4 km s" 1 Mpc" 1 , 
n m = 0.262, fU = 0.738 (Komatsu et al. 2010). 



APPENDIX B: ALGORITHM OF 
CALCULATING 5d FOR MODEL DISCUSSED 
IN SEC. 4.2 

The algorithm of calculating 3d is as follows 

(i) First the radius of a void R v is generated from the 
Gaussian distribution with the mean of 12 Mpc and the 
standard deviation of 2 Mpc. 

(ii) Density within the void Q v is generated from the 
Gaussian distribution with the mean 0.2fi m and a — 
0.27f2 m . If a generated in this way Q v is lower than 0.01f2 m 
then the generation is repeated. If after 6 times it is still less 
than 0.01Sl m then Q v is generated for a uniform distribution 
between and 0.01f2 m . If Q v Q m then Q v is generated one 
more time. If after 6 times £l v ^ Q m then its value is chosen 
for a uniform distribution from 0.85f2 m to Q m . 

(iii) Density of the surrounding shell Q s is generated from 
the Gaussian distribution with the mean of 1.75Sl m and a — 
0.7f2 m . If Q s ^ flm then its value is generated again. If after 
6 times il a ^ Q m then its value is generated for a uniform 
distribution between 1.75f2 m and 1.95f2 m . 

(iv) The condition, that the structure is compensated im- 
plies that the radius of the whole structure is 

(v) The angle at which the light ray enters the structure 
is generated using 4 different methods - for details see Sec. 

(vi) The evolution of 8 (at a fixed point) was calculated 
using (|ATj) . 

(vii) Integral ([SJ) is solved from Xi to Xf (where Xi is the 
comoving coordinate of the entry point and Xf the point 
where the light ray exits the structure). 

(viii) Steps (i)-(vii) are repeated so © is solved from 
X = to Xe- 

If instead of the Gaussian, the uniform PDF is used 
(steps (i)-(iii)) then the PDF of 5 is not log-normal PDF as 
presented in Fig. [3] 

The cosmological parameters are the same as derived 
from the 7-year WMAP data: H = 71.4 km s" 1 Mpc" 1 , 
Q m = 0.262, Q A = 0.738 (Komatsu et al. 2010). 

This paper has been typeset from a Tf^X/ DTJ5X file prepared 
by the author. 
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